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Abstract. We discuss the existence of a mass-metallicity relation among galactic halo globular clusters. The 
lack of any luminosity- metallicity correlation in globular cluster systems has been used as an argument against 
self-enrichment models of cluster formation. We show that such a relation is statistically present among the 
galactic Old Halo globulars. This observational correlation implies that the least massive old clusters are the most 
metal-rich. This is in contradiction with the idea that, if globular clusters were self-enriched systems, the most 
metal-rich clusters would also be the most massive ones. We further show that this ant i- correlation is as predicted 
by self-enrichment models. 
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1. Introduction 

Galactic Globular Clusters are among the very first 
bound structures that formed in what later became the 
Galaxy. As such, individual globular clusters and the 
Globular Cluster System provide us with important 
insights concerning the age, the formation and the early 
evolution of the Galaxy. 

This paper is the third of a series devoted to the study 
of a formation scenario of galactic halo globular clusters, 
namely the self-enrichment hypothesis, which develops 
the Fall & Rees ( 1985| ) cluster formation model. The 



model is detailed in Parmentier et al. ( 1999 ) (hereafter 
Paper I) and a summary is provided in Parmentier et 
al. (|2000D (hereafter Paper II), Sect. 2. 



The model assumes that primordial cold clouds embedded 
in a hot and diffuse protogalactic background (Fall & 
Rees 1985) are the gaseous progenitors of galactic halo 
globular clusters, that is, this model assumes baryon 
assembly predates star formation. Our model explores 
the ability of these proto-globular cluster clouds to retain 
the ejecta of a first generation of zero-metal abundance 
stars, born in the central regions of the clouds. When 
the massive stars explode as Type II supernovae, they 
chemically enrich the surrounding pristine interstellar 
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medium and trigger the expansion of a supershell in 
which a second generation of nonzero-metal abundance 
stars may form. The aim of a self-enrichment scenario 
is therefore to explain both the formation of a globular 
cluster and the origin of its metal content. 

One of the key parameters of this class of model is 
the external pressure exerted by the hot protogalactic 
background on the proto-clusters. The higher the pres- 
sure is (i.e. the deeper the proto-cluster is located in 
the protoGalaxy in the simplest implementation of the 
model), the smaller its mass is, the higher its metallicity 
will be (see Table 1 of Paper I) . 

An in-depth discussion of the ensuing Galactic metallicity 
gradient is presented in Paper II. We show that, when 
combined with a pressure profile scaling as Ph oc D~ 2 , 
where Ph is the hot protogalactic background pressure 
and D is the galactocentric distance, the model is con- 
sistent with the metallicity gradient observed for the Old 
Halo globular cluster system. 

There are three aspects of globular cluster forma- 
tion which self-enrichment models must specifically ad- 
dress. The disruptive effects of supernovae, and the inter- 
nal chemical homogeneity are discussed in Paper I. This 
paper considers the third aspect, the extent to which a 
mass(luminosity) -metallicity relation is expected and ob- 
served. The Galactic globular cluster system is used as a 
specific example. 
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Fig. 1. Mass- luminosity relati on fo r the 56 globular clus- 
ters of the Pryor & Meylan ( 1995 ) compilation and the 
corresponding least-squares fit 



2. Self-enrichment and a mass-metallicity relation 

The lack of any obvious correlation, in any globular clus- 
ter system, between the mass (or the luminosity) and the 
metallicity of individual globulars is often used as an argu- 
ment against the self-enrichment hypothesis. Indeed, were 
one to assume that gravitational potential gradients dom- 
inated mass loss, the most massive objects would be bet- 
ter able to retain their metal-enriched supernova ejecta, so 
that metal abundance should increase with cluster mass in 
case of self-enrichment. Before adressing the discussion of 
a luminosity-metallicity relation, we would like to dismiss 
this idea that more massive clusters would be more metal- 
rich in the case of self-enrichment. If a more massive ob- 
ject is indeed better able to retain more supernova ejecta, 
this larger amount of metallic ejecta is mixed with a larger 
amount of gas. Therefore, no firm conclusion can be drawn 
concerning the resulting metal abundance (or metallicity), 
i.e. the ratio of the two increased quantities. It is the 
fractional efficiency of gas retention which is important. 
Most importantly, though, mass loss in this class of mod- 
els is determined by external gas pressure and not by the 
pressure equivalent of the gravitational potential gradient. 
This means that the absence of a mass-metallicity relation, 
in the sense that the most massive globul ars w ould also 
be th e mos t metal-rich (e.g. McLaughlin 1997 , Barmby 
et al. 



200Q) , can not be considered as evidence against 



the self-enrichment hypothesis. In sharp opposition with 
these statements, the self-enrichment model we develop 
foresees a mass-metallicity relationship in the sense that 
the most metal-rich proto-globular clusters are the least 
massive ones. 

Unlike globular clusters, dwarf galaxies exhibit well- 
defined correl ation s betwe en lum inosity and metallicity 
(e.g. Gilmore |2000| , Mateo |2000| ) such that the dimmest 
ones are the most metal weak. The standard explana- 
tion for this correlation being sclf-cnrichmcnt in the pres- 
ence of galactic winds which are limited by gravitational 
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Fig. 2. Mctallicity-luminosity diagram for the whole 
galactic globular cluster sytem. The dashed line at 
[Fe/H]=— 0.8 represents the generally assumed metallic- 
ity limit between the halo and the bulge/disk subsystems. 
The different types of globulars are marked by different 
symbols. Old Halo and Younger Halo clusters are respec- 
tively represented by open and full circles. The crosses la- 
bel the Sgr globulars. The open squares ("noHBR" group) 
stand for the halo globular clusters for which the horizon- 
tal brach morphology index is not given in Harris (1996). 
The full triangle and the asterisks respectively represent 
u) Cen and the bulge/disk clusters 



potential grad ients (Dekel & Silk [1986|), Djorgovski & 
Meylan ( 1994 ) conclude that globular clusters c annot 
be self-enriched systems. However, Dekel & Silk ( 1986 ) 
point out that the dwarf galaxy observed luminosity- 
metallicity relation can be successfully explained only if 
the gaseous proto-galaxies are embedded within domi- 
nant halos of dark matter. While there is indeed clear 
evidence of the presence of such halos around dwarf galax- 
ies (Mateo 1996), this is not the case for globular clus- 
ters (Moore [1996] Meylan & Heggie |1997| ). Therefore, the 
Dekel & Silk model, built for dwarf galaxies, can cer- 
tainly be not extrapolated to globular clusters. Moreover, 
dwarf galaxies and GCs have undergone very different star 
formation histories: their respective star formation rate 
and duration differ by, at least, an order of magnitude 
(Gilmore 2000). Dwarf galaxies also exhibit metallicity 
spreads, often larger than ldex (Mateo 2000), in marked 
contrast with the chemical homogeneity of globular clus- 
ters. Considering these many differences, the comparison 
between globular clusters and dwarf galaxies therefore ap- 
pears irrelevant. 

In searching for a luminosity-metallicity relation in 
the Galactic globular cluster system, it should be kept 
in mind that, while the observed quantity is the luminos- 
ity, the physical quantity of interest is the mass. Figure [l] 
represents the relation between the mass and the abso- 
lute visual magnit ude f or the 56 globular clusters of the 
Pryor & Meylan ( 1993 ) mass compilation. The M v val- 
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ues come from the McMaster Catalogue (Harris 1996, up- 
dated 1999). The scatter superimposed on the correlation, 
of the order of ctm v —0.6, is equivalently the variations of 
the mass-to-light ratio from cluster to cluster. This ranges 
from ~1 to ~4 (Pryor & Meylan 1993) and reflects pos- 
sible differences in the initial mass function and the dy- 
namical evolution of the clusters n. Therefore, any mass- 
metallicity correlation will be, at least partly, washed out 
in the corresponding M„-[Fe/H] plot. This effect is illus- 
trated in Fig. ^, the metallicity-luminosity diagram for the 
whole globular cluster system (the [Fe/H] values are taken 
from the McMaster Catalogue). Also plotted are t he cor - 
responding error bars in [Fe/H], ± 0.15dex (King |l999|) , 
and in M v , ± 0.6 dex from Fig. [y, if the latter is considered 
to be a mass indicator. The size of the M v errorbars (re- 
flecting the different luminosities that a globular cluster 
with a given mass but varying mass-to-light ratios may 
exhibit) is clearly not negligeable compared to the size 
of the observed distribution, the dispersion of the best- 
fitting gaussian to the galactic globular cluster luminosity 
function being ~1.2 (Harris 1991 ). 

Unfortunately, determination of the physical quantity 
of interest, i.e. the relative masses of the globular clus- 
ters at their fo rmati on, is still uncertain at least by a fac- 
tor 2 (Meylan 2000 ). For instance, the use of single-mass 
King models is a simplification which tends to underes- 
timate cluster mass (Ashman & Zepf 1998 Mcrmilliod 
200C). Therefore, one of the key points in the search for 



a mass-metallicity correlation is to use an homogeneous 
set of globular cluster masses in order to limit additional 
scatter in the (logMgc [Fe/H]) plot. We use the globular 
cluster mass compilation computed by Pryor & Meylan 
( 1993 ): this compilation is the most complete set of glob- 
ular cluster masses computed with an internally consistent 
family of multi-component King-Michie models. 

Another source of scatter in the luminosity(mass)- 
metallicity plot is introduced by the various origins of the 
Galactic globulars. Indeed, evidence has now accumulated 
that the Galactic globular cluster system does not consist 
of globular clusters with a single origin. While the ma- 
jority of globular clusters in the halo are ol d, wit h a re- 
markably small age spread (Rosenberg et al. 1999), there 



is a small subset, particularly among the more metal-rich 
clusters, with inferred ages of several Gyr younger than 
the dominant old population. These younger globular clus- 
ters are either clusters being/having been accreted by the 
Galaxy recently or metal-rich clusters associated with the 
bulge/disk subsystem. These clusters being significantly 
younger, their formation is not expected to be taken into 
account by our self-enrichment model, which deals with 
globular clusters whose gaseous progenitors have a pri- 
mordial composition. 
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Fig. 3. [Fc/H] vs log Mac diagr am fo r the globular clus- 
ters studied by Rosenberg et al. ( 1999 ). The old (OC) and 
the younger (YC) globular clusters exhibit different distri- 
butions in the diagram relatively to the boundary foreseen 
by the self-enrichment model. Also shown is the tidal de- 
struction zone for the less massive halo globular clusters 
located inside the solar circle 

confirms the globular cluster system subdivisions early 
suggested by Zinn ( 1985| , 1993 ). From the point of view 
of the metallicity distribution, the Galactic globular clus- 
ter system is composed of two subpopulations, a metal- 
poor halo group and a metal-rich, centrally concentrated, 
bulge (or disk) group (Zinn |l985| ). Further more, the halo 
subsystem itself includes an Old Halo, made of globular 
clusters perhaps born in situ, during the collapse of the 
protogalactic cloud, and a Younger Halo likely made of 
globulars later stripped from neighbouring dwarf galaxies 
(see Paper II for a review of these evidences). 

Since our self-enrichment model deals with globular 
clusters whose progenitors were embedded in the hot 
phase of the protogalactic cloud and whose gaseous mate- 
rial was pristine, it is not expected to apply to the Younger 
Halo group, the presumed accreted component of the halo, 
nor to the bulge clusters. Thus, in what follows, we fo- 
cus eithe r on the coeval and old sample of Rosenberg et 
al. (1999) or on the Old Halo defined by Zinn. 



3. Comparison of the model with the observations 

Our self-enrichment model suggests the existence of an 
anti-correlation between the mass M of a proto-cluster 
and the metallicity [Fe/H] reached at the end of the self- 
enrichment process, in the sense that the least massive 
proto-clusters create the most metal-rich globular clusters 
(see Table 1 of Paper I): 



The age spread highlighted by Rosenberg et al. ( |1999[ ) also [Fe/H] = 4.3 — log M 



(1) 



1 Part of the scatter also originates in the uncertainties of 
globular cluster mass computation, since it relies on model as- 
sumptions such as the isotropy of velocity distibution and a 
maxwellian distribution of velocities. 



However, this (anti-) correlation applies to the gaseous 
progenitors; the mass-metallicity relation observed among 
the studied sample of clusters, if any, should only be 
a relic of Eq. (fil). For instance, the —1 slope will be 
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Fig. 4. [Fe/H] vs log Mqc plot including Old Halo and 
Younger Halo subgroups (49 globular clusters). 



conserved only if there is a universal and constant star 
formation efficiency for the second stellar generation 
which forms the majority of the stars from the chem- 
ically enriched gas swept up from the first generation 
supernovae. Since the mass M of a gaseous progenitor is 
an upper limit for the mass Mac of the globular cluster 
formed, Eq. (|l|) delimits a permitted area in the (logMac, 
[Fe/H]) plot: all the data should be located to the left of 
Eq. (1) (plain curve in Fig. || and ||). Figure |] represents 
Eq. (1) together with the globular clusters for which the 
age an d the mass are respectively prov ided in Rosenberg 
et al. fll999| ) and Pryor & Meylan ( |l993[ ). Obviously, 
the two Rosenberg et al. ( 1999 ) groups (old clusters: 
open symbols; younger clusters: filled symbols) behave in 
a different way compared to our self-enrichment mass- 
metallicity relation. While the old, coeval and metal-poor 
GCs are all located in the permitted area of the plot, 
i.e. their mass-metallicity diagram is consistent with the 
self-enrichment of primordial gaseous progenitors, half of 
the young clusters, either presumed accreted or belonging 
to the bulge subsystem, are located in the forbidden area 
of the diagram, i.e. as expected, their formation cannot 
be accounted for by the self-enrichment model. In Fig. ||, 
we also represent three of the most massive globulars 
(filled triangles), namely w Cen, M 54 and NGC 5824. 
Their location in the forbidden part of the plot points to 
a different star formation history. This is not surprising 
since, at least in the case of u> Cen and M 54, an intrinsic 
abundance spread is seen. M 54 is of course a member of 
the Sagittarius dwarf spheroidal galaxy, and is not (yet) 
a Galactic globular cluster. 



Figure || is also clearly depleted in low-mass globular 
clusters (log Mqc < 4.8). However, at a galactocentric 
distance smaller than 8kpc, these low-mass clusters are 
not expected to survive more than a Hubble time (see the 
"survival triangle" in the mass vs half-m ass radius dia- 
gram defined by Gnedin & Ostriker 1996 , their Fig. 20a). 
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Fig. 5. [Fe/H] vs log Mac plot for the Old Halo subgroup 
(37 globular clusters). The self-enrichment model (plain 
curve) defines a permitted area (left part of the plot) in 
which most of the observational points are located. The 
correlation between the globular cluster masses and their 
metallicities is particularly striking for [Fe/H]> —1.8 



The vast majority of the globular clusters located at 
these galactocentric distances, i.e. log D< 8 kpc, exhibit a 
metallicity higher than [Fe/H] =—2. The depletion zone, 
represented by the box in Fig. || and ||, is therefore not 
surprising and corresponds to the tidal destruction of 
these low-mass clusters. The globular clusters used in our 
Paper are therefore no more than a surviving sample. 
The distance a given cluster lies to lower masses from 
the model upper bound is, to first order, a measure 
of the star formation efficiency of cluster formation. A 
"typical" surviving cluster lies a factor of order 5 below 
the bound, suggesting an efficiency factor of order 20%. 
As noted above, however, lower mass clusters will have 
preferentially failed to survive until today, so that this 
value is an upper limit. Star formation efficiencies in 
the range from a few to a few tens of percent seem 
appropriate for most clusters. Only the few percent of 
clusters which are the most massive require star formation 
efficiencies in excess of unity, and so are inconsistent with 
this formation model. Interestingly, these very massive 
clusters are those which show internal abundance spreads, 
which are themselves direct evidence for self-enrichment 
during cluster formation. 

In order to increase our sample and to look for a 
surviving correlation between the mass and the metallic- 
ity, we also consider the Old Halo subgroup (Zinn |1993 ). 
As for the metallicity gradient (see Paper II), an Old 
Halo/ Younger Halo separation is fruitful. Figure [I] shows 
a plot of [Fe/H] versus mass for the 49 halo globular clus- 
t ers w hose mass has been computed by Pryor & Meylan 
( |1993|) : there is no correlation between the mass and the 
metallicity, the linear Pearson correlation coefficient being 
-0.15. 
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Considering the Old Halo group only (Fig. ||), as stated 
in the previous Section, a weak correlation between the 
logarithm of the mass of the globular clusters and their 
metallicity emerges. The linear Pearson correlation coeffi- 
cient improves to a value of —0.35, with a corresponding 
probability of correlation of 96.92%. Moreover, most of 
the Old Halo globulars are located in the permitted area 
of the plot. | 

4. Conclusions 

We present the mass-metallicity relation foreseen by our 
self-enrichment model. At first glance, such a globular 
cluster formation scenario is disproven by the lack of any 
obvious correlation between the luminosity and the metal- 
licity of globular clusters, a s claim ed by previous a uthor s 
(e.g. Djorgovski & Meylan 



1994 



Ashman & Zepf [L99 



However, we stress here that there are numerous sources 
of scatter between the theoretical (mass, [Fe/H]) relation, 
applying to the gaseous progenitors of globular clusters, 
and the observed (luminosity, [Fe/H]) plot, applying to 
globular clusters. These scatter sources are, for instance, 
the star formation efficiency with which the globular clus- 
ter stars form out of the proto-globular cluster cloud, 
the mass losses undergone by globular clusters with time 
(see the tidal tails exhibited by some clusters), the varia- 
tions in the mass-to-light ratio from one globular cluster 
to another. We also caution that the search for a mass- 
metallicity correlation should be restricted to a given glob- 
ular cluster subpopulation, namely the Old Halo group. 
Despite the numerous sources of scatter, the globulars of 
this group are characterized by a boundary in the mass- 
metallicity diagram and by a correlation in the sense ex- 
pected by simple pressure-bounded self-enrichment mod- 
els. Self-enrichment models remain a viable hypothesis for 
galactic halo globular cluster formation. 
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